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Abstract 

Mixed dark matter models with one low-mass (e.g. 2.4eV) neutrino flavor are 
shown to be in good agreement with observations if the neutrinos have non¬ 
zero chemical potentials. This agreement holds (except for the problem with 
a low-age Universe) even for high values of the Hubble-parameter. Massless 
neutrinos with non-zero chemical potentials may reconcile cold dark matter 
with observations. Some fine-tuning is required to avoid problems with Big 

Bang nucleosynthesis. 
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The mixed dark matter (MDM) model has been very successful in reproducing the large 
scale structure observed in the Universe [1], and in fact is among the only models that can 
at the same time fit both the very large scales, as measured by COBE [2], and the studies 
of galaxy correlations, as recently summarized by Peacock and Dodds [3]. Very nice fits 
[4] are obtained if one neutrino of mass near 4.8eV (or even better, 2 neutrinos of mass 
near 2.4eV) is responsible for roughly 20% of the density in the Universe, with another 80% 
(apart from a minor fraction of baryons) in terms of an unknown form of cold dark matter 
(CDM). The relevance of such models is significantly strengthened if preliminary indications 
from oscillation experiments [5] of a 2.4eV muon neutrino mass are confirmed. 

However, standard MDM models only fit the data if the expansion rate of the Uni¬ 
verse has a definite value. In terms of the Hubble parameter, Hq, this value is Hq 
50 kms“^ Mpc“^, or ho = ifo/100 kms“^ Mpc“^ 0.5. Several recent investigations [0], 
while still marred by possible large systematic uncertainties, tend to give much higher val¬ 
ues, in the range of ho ~ 0.7-0.8, which would be inconsistent with the MDM-models. (Such 
a high Hubble parameter would also make the Universe suspiciously young unless the cosmo¬ 
logical constant is non-zero. An U = 1 Universe has an age of only to = 6.5/io^Gyr, giving 
to < 9.3Gyr for ho > 0.7. Standard MDM-models are however in conflict with large scale 
structure data (regardless of the choice of m,^) even for ho = 0.6, where the age problem is 
not severe. In the present investigation we shall not dwell further on the age problem, but 
concentrate on reproducing the observed large scale structure). 

The inconsistency between MDM and ho > 0.5 is related to the coupling between neutrino 
mass, neutrino density contribution, and neutrino free-streaming. As pointed out in the 
present investigation, mixed dark matter models can be reconciled with ho > 0.5 if this 
coupling changes, as in the case of a non-zero neutrino chemical potential [7]. Due to the 
preliminary experimental indications of a 2.4eV we shall focus in particular on a scenario 
where one such flavor is the hot component, but similar conclusions can be drawn for other 
choices of mass and number of flavors. 

Neutrinos decoupling while ultra-relativistic are described by the momentum distribution 
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with T = Tuan I a, fJ- = iido-dIo., where T, a, and fi denote temperature, scale factor, and 
chemical potential, and subscript D indicates the value at decoupling. 

In the detailed investigations described below we have integrated the full neutrino dis¬ 
tribution function to get the density, free-streaming etc. Some approximate expressions in 
the ultra-relativistic and non-relativistic limits are useful, however. 

Assuming neutrino-antineutrino equilibrium (/x = —/x) the present neutrino contribution 
(including u) to the cosmic density in units of the critical density is 

= 2.382 X ^77X2.4 [ 77 ^ + vr^77 + 6A2 (—t?)] • (2) 



In this relation 77 = /x/T and 77 x 2.4 = 77x,^/2.4 eV. The neutrino is assumed to be single- 
handed. The function F 2 is 

00 p-ir + i)v 

r,(-,) = r(3)i:(-ir^. (3) 

taking the value 3r(3)^(3)/4 for 77 = 0. Thus for non-degenerate neutrinos, = 2.58 X 
10“^/xq ^ 77 x 2 . 4 , giving 0.21 for a 4.8eV neutrino or two 2.4eV neutrinos and ho = 0.5. 

In the ultra-relativistic regime a degenerate neutrino-antineutrino flavor contributes to 
the density like one non-degenerate flavor plus an additional 
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flavors. A problem with such a model is therefore the increase in relativistic density during 
Big Bang nucleosynthesis, leading to an overproduction of helium. To avoid this one needs 
for instance a slight degeneracy in electron neutrinos to push the weak interactions more in 
the direction of protons. Several authors have demonstrated that nucleosynthesis data only 
constrain the ratio that the observations are in agreement with a large 

and/or chemical potential (see [8] for a review). In contrast to the cosmic net baryon 
number, which is known to be small, there are no strong observational or theoretical con¬ 
straints on the lepton number, which is a sum of contributions from electrons and 3 neutrino 
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flavors, where only the electron term is known to be small (from charge neutrality). The 
total lepton number could be large or it could be comparable to the baryon number. If it is 
comparable to the baryon number, this could occur because all neutrino terms were small, 
but it could also be due to cancellations. Even if a (GUT scale?) neutrino chemical poten¬ 
tial is small, non-zero effective chemical potentials could appear later from particle decays 
[9]. While non-zero chemical potentials may seem like an ad hoc solution to the structure 
formation problem, it is not in conflict with any fundamental principles or observations. 

We note in passing that a massless (or very low-mass) degenerate neutrino could be 
responsible for the increase in relativistic degrees of freedom needed in the scenario of Do- 
delson, Gyuk and Turner [10] to reconcile a pure GDM scenario with observations. 

Finally it is worth noting that the mean momentum of degenerate neutrinos is higher 
than that of non-degenerate neutrinos (3/x/4 in the limit rj ^ oo versus S.lblT for rj = 0). 

The relations above show the qualitative effects of neutrino degeneracy to be the follow¬ 
ing: A larger density contribution for fixed m,^, allowing a smaller neutrino mass to explain 
a given HDM-fraction for fixed ho, or permitting a larger Hubble parameter for fixed m,^. At 
the same time, the free-streaming for fixed rrii, and hg increases somewhat. Detailed calcu¬ 
lations summarized below show, that the final density perturbation spectrum is consistent 
with observations for a scenario with one degenerate 2.4eV neutrino (for instance the muon 
neutrino, assuming then an unstable and any choice of Hubble parameter in the interval 
0.5 < ho < 1. This in contrast to non-degenerate neutrinos, where ho > 0.5 is inconsistent 
with observations for any choice of neutrino mass. 

We calculate the evolution of an initial density perturbation in a flat (H = 1) Universe 
using linear theory and follow the approach by Bond and Szalay [11]. We assume one 
massive neutrino species with mass rrii, = 2.4 eV and a non-zero chemical potential as the 
hot dark matter (HDM) component, leaving two massless neutrino species. The fraction 
of baryons in the Universe is only a few per cent, so to simplify calculations we neglect 
them completely when calculating the evolution of density perturbations. As our aim is 
to demonstrate the general features of a degenerate neutrino model rather than present 
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detailed results for specific parameters the errors introduced thereby are not crucial. The 
CDM fraction is therefore just flcDM = 1 ~ 

We include the gravitational effect of photons in the tight-coupling approximation [11], 
neglecting photons after recombination. Only adiabatic density perturbations are consid¬ 
ered, with initial conditions [11] hcDM = The scale factor a is normalized such 

that a = 1 when the massive neutrinos become non-relativistic. We integrate the differential 
equations given in [11] (properly modified to allow degeneracy) from a = 0.005 and until 
o- = Unow — 4000. The scale factor today Unow depends on the chemical potential since the 
mean momentum of the neutrinos grows with r]. We assume 77 > 0 so that antineutrinos 
have fj = —7] < 0. This means that the neutrinos will have a larger mean momentum 
than in the usual rj = 0 case, and therefore more free-streaming so that more of the initial 
density perturbations will be erased, whereas the antineutrinos will free-stream less than 
usual. It is therefore practical to evolve the massive neutrino perturbations as two distinct 
perturbations 61, and 6p. For r) > 2 neutrinos greatly outnumber antineutrinos. 

The very large scale density perturbations, A > SOO/ig ^ Mpc enter the horizon well after 
the massive neutrinos have become non-relativistic. Thus the perturbations on these large 
scales are not damped by any free-streaming and just grow like 6 oa when the Universe is 
radiation dominated, and as h oc a when the Universe shifts to be dominated by matter. The 
small scale perturbations enter the horizon when the massive neutrinos are still relativistic 
and therefore they are damped by free-streaming. The CDM perturbations grow at a lower 
rate than hcDM oc a as long as the HDM perturbations 61 , and 6 p are smaller than hcDM! the 
more homogeneous HDM distribution suppresses the growth of CDM. 

We calculate the transfer function T[k) for ten wavenumbers ranging from k = 4: X 
10“^/ioMpc“^ to k = 2/ioMpc“^ (the corresponding wavelengths are A = 27r/A;; all wave¬ 
lengths and wavenumbers are comoving with values given as they would appear to be to¬ 
day). The initial power spectrum is assumed to be of Harrison-Zel’dovich type as predicted 
by many theories of inflation. The power spectrum is given as P{k) = AkT^{k), where 
the normalization factor A is determined by the quadrupole anisotropy of the microwave 
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background radiation measured by COBE [2], For this we have taken the value Q = 17/xif. 
In order to compare directly with the large scale structure data extrapolated to the linear 
regime by Peacock and Dodds [3] we follow their notation and calculate the dimensionless 
power spectrum, 

^ ( 5 ) 

which we henceforth will refer to as the power spectrum. A^(A;) can be described as the 
contribution to the fractional density variance per logarithmic interval in k. 

With the chemical potential as a free parameter, we calculate the power spectrum for a 
range of values for hg from 0.5 to 1.0 adjusting the value of rj for each value of hg in order 
to obtain a reasonable fit to the large scale structure data points [3]. In Figure 1 we show 
the calculated power spectrum for six models with values of hg ranging from 0.5 to 1.0; the 
corresponding values of 77 , and shape parameter P (to be defined later) are shown in 
Table 1. 

The power spectrum A^(A;) is plotted as a function of the wavenumber in units of the 
Hubble parameter hg. Almost the same shape of the power spectrum is obtained for the six 
different values of hg. 

If the density fluctuations are Gaussian (those generated by inflation are), all observable 
quantities of large scale structure can be directly calculated from the power spectrum as long 
as the perturbations are small, that is for only linear evolution. Thus if two power spectra 
are equal they would predict the same amount of large scale structure today. Within a few 
per cent the best-fitting models have rj = 4.7 hg — 0.53. (We note again, that even though we 
And very good agreement with the data for 0.5 < ho < 1, the age problem probably forces 
hg < 0.7). 

Figure 1 also shows the standard CDM model and the MDM power spectrum for one 
2.4eV neutrino flavor with rj = 0 and hg = 0.5. Both models are clearly inconsistent with 
the data. Increasing hg would make the curves bend the wrong way, making the fit even 
worse. 
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In order to further test if the calculated models can fit the properties of the real Universe, 
we derive some observable quantities from the power spectra. The rms mass fluctuation in 
a top-hat sphere of radius r = S/ig^Mpc, CS) is found to be around 0.75 in accordance with 
other MDM models. From observations of galaxy surveys the rms fluctuation in galaxy 
number density is usually found to be Us^gai = 1 , so galaxies are more clustered than mass, 
corresponding to a bias parameter of 6 = crs^gai/crs = 1-3. 

Another somewhat more direct measure of the mass distribution is the peculiar bulk 
velocity in spheres, which for a radius of SO/ig^Mpc is observed to be V 50 = 335 ± 80km/s 
[4]. We And a model-value between 370km/s and 390km/s for V 50 in fine agreement with 
observations. 

The number density of clusters of galaxies, fVciustj with a mass M > can be 

found using the Press-Schecter approximation [12]. We follow [4] and [13] using a Gaussian 
Alter and the density contrast needed for collapse assumed to be 6 c = 1.5. Observations 
give fVciust = 4.0 ± 2.0 x 10 “^/iqMpc“^ [4], whereas we get number densities from 6.4 to 
10.4 X 10 ^/iqMpc ^ for the power spectra shown in Figure 1 when ho ranges from 1.0 
to 0.5. Thus our models slightly overproduce clusters, particularly for hg = 0.5. Proper 
inclusion of baryons will reduce the power on cluster scales and below because initial density 
perturbations in baryons are not allowed to grow before recombination. This would reduce 
Afciust- It is also possible to reduce the power on large scales by choosing a lower normalization 
(there could be some gravitational wave contribution to the COBE anisotropy), and at the 
same time increase the power on small scales by choosing a smaller value of 77 . 

The scale where the auto-correlation function crosses zero is observed to be larger than 
40/io^Mpc [4]. We And a consistent value of 56/io^Mpc with very little variation. In general 
MDM models give a better At than standard CDM, which predicts zero crossing already at 
36/io^Mpc [4]. 

Thus MDM models with degenerate neutrinos At the observed large scale structure rather 
well. On smaller scales of A ~ l/i“^Mpc or less non-linear eflects are crucial, and comparison 
with observations less reliable. One comparison which has been much discussed recently is 
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with the observed number density of damped Lyman-alpha systems (DLAS), which constrain 
the mass fraction of collapsed baryons at high redshifts [14], MDM models form the structure 
quite late and have some difficulty in forming enough collapsed structures to explain the 
DLAS [13,15]. From the observations which are quite uncertain at high redshifts, the mass 
fraction of collapsed baryons at 2 ; = 3 — 3.5 is determined to be flgas = 6.0 ±2.0 X 10“^. From 
the models ho = 0.5 and ho = 1.0 (the ones which differ most, particularly at smaller scales) 
we find values of flgas to be at least 3.9 X 10“^ and 6.5 X 10“^ respectively, by following the 
methods described in [13]. (For the large wave numbers needed to calculate flgas, we have 
approximated the evolution of the power spectrum by the growth of the CDM component 
in a flat Universe with a uniform background density of neutrinos The mass fraction 

of collapsed baryons was obtained assuming Ufc = 0.075 and assuming that all the gas is 
neutral, but some of the gas could be ionized or even removed by early star formation, thus 
the observed value of Ugas should here be taken as a lower limit. As we have not included 
baryons in the calculation of the power spectrum, the limits on Ugas should be taken only 
as an indication, that it is possible to make our MDM models with a chemical potential 
consistent also with the constraints from DLAS. 

What is the dominant effect that modifies the initial power spectrum? The power spec¬ 
trum is suppressed on small scales, because density perturbations inside the horizon can 
not grow before the Universe becomes matter dominated at time teq- Before this time all 
the neutrinos are relativistic (the massive neutrinos become non-relativistic just around teq)- 
The shape of the power spectrum at tgq is roughly determined by the scale that crosses the 
horizon at that time, Agq. On scales larger than Agq, the transfer function T{k) = 1, whereas 
on scales smaller than Agq the transfer function decreases roughly as A^ for pure CDM, 
and even stronger in MDM when neutrino perturbations are erased by free-streaming. The 
scale Aeq can be written as Aeq ~ lOMpc(Uoho) ^(5'*/3.36)^/^ [ 10 ], where 5 * is the effective 
number of particle degrees of freedom. The shape of the power spectrum is determined by 
(hg ^Mpc/Aeq), so one can define a shape parameter, F = Uoho(5'*/3.36)“^/^ [10]. In standard 
MDM models one has Ug = 1 and ho = 0.5 with three neutrino flavors giving F = 0.5. Stan- 
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dard CDM can be made to fit the large scale structures observed by modifying the shape 
parameter F, the required value of F 0.25 [3] can be obtained by adding more relativistic 
degrees of freedom increasing 5 * [ 10 ] (as noted earlier, a massless, degenerate neutrino might 
do exactly that). 

In the MDM models shown in Figure 1, the increase of 77 needed to fit the observed 
power spectrum for high values of ho increases the effective number of neutrino species in 
the relativistic era by given in Equation (4). This again increases 5 * = 3.36 + 0.454AfV,^ 
and makes the shape parameter F stay almost constant at around 0.5 as needed for MDM 
models with ~ 0.25; see Table 1. 

The computed models are of course complex, and can only roughly be characterised by 
the shape parameter, in particular because neutrinos become non-relativistic very close to 
feq- The free-streaming properties also differ a bit. When 77 is high the number density of 
massive neutrinos is higher although stays roughly the same, and the mean momentum 
increases a little. Still, the dominant effect that modifies the power spectra, is the increase 
in radiation energy density delaying feqj which cancels out the increase in Hubble parameter. 

We have assumed only one massive neutrino species, being the muon neutrino with a 
mass of 2.4 eV. There are of course other possibilities. Primack et al. [4] favor their Ct/^DM 
model, which has 2 massive neutrino species (muon and tan) both with masses of 2.4 eV. 
That particular model makes a remarkably good fit to all the observed values of both large 
and small scale structures if hg ^ 0.5. It would be equally easy to make MDM models with 
two massive neutrinos and a medium/high value of the Hubble parameter fit the observed 
large scale structure by introducing a chemical potential for one or both of the neutrino 
species. Similarly for other values of the neutrino masses. Also, we have only considered 
a scenario where the massive was degenerate. The extra relativistic degrees of freedom 
could have been produced by a degenerate, massless neutrino as well. 

To summarize, we have demonstrated how neutrino degeneracy makes MDM models vi¬ 
able for one low-mass neutrino flavor, even for medium/large values of the Hubble parameter. 
It would be necessary to perform more precise calculations of the power spectrum in order 
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to determine exactly how well MDM models with a chemical potential can be made to fit all 
the observational constraints. No doubt this will be done if the neutrino mass detections are 


confirmed and the evidence for a large value of the Hubble parameter continues to increase. 
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FIGURES 

FIG. 1. The dimensionless power spectrum A^(A:) plotted as a function of wavenumber in units 
of /ioMpc“^. The MDM power spectra for six different choices of Hubble parameter ranging from 
0.5 to 1.0 are shown, all normalized to the COBE quadrupole anisotropy. For each power spectrum 
we have chosen the value of rj such that a good fit to the observed linear power spectrum is obtained; 
the data points are taken from [3]. The corresponding values of rj are given in Table 1. Also shown 
for comparison is the standard CDM power spectrum, and the MDM power spectrum with ho = 0.5 
and rj = 0 for one 2.4eV neutrino flavor (solid curve close to CDM). 
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TABLES 

TABLE 1. Corresponding values of Hubble parameter, degeneracy, neutrino density fraction, 
and shape parameter, for the MDM models of Figure 1. 


ho 

V 

^1/1/ 

r 

0.5 

1.75 

0.235 

0.45 

0.6 

2.3 

0.239 

0.50 

0.7 

2.8 

0.245 

0.54 

0.8 

3.3 

0.257 

0.56 

0.9 

3.7 

0.257 

0.58 

1.0 

4.1 

0.261 

0.59 
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